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ELEMENTS OP THE WING SECTION THEORY AND OF THE WING 
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By Max M. Munk:. 



SUMMARY. 



The following paper, prepared for the National Advisory Committee for Aeronautics, 
contains those results of the theory of wings and of wing sections which are of i mm ediate 
practical value. They are proven and demonstrated by the use of the simple conceptions of 
"kinetic energy" and "momentum" only, familiar to every engineer; and not by introducing 
"isogonal transformations" and "vortices," which latter mathematical methods are not 
essential to the theory and better are used only in papers intended for mathematicians and 
special experts. 
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I. THE COMPLEX POTENTL4L FUNCTION. 

1. I have shown in the paper, reference 1, how each air flow, considered as a whole, 
possesses as characteristic quantities a kinetic energy and a momentum necessary to create it. 
Many technically important flows can be created by a distribution of pressure and they then 
have a "velocity potential" which equals this pressure distribution divided by the density 
of the fluid with the sign reversed. It is further explained in the paper referred to how the 
superposition of several "potential flows" gives a potential flow again. 

The characteristic differential equation for the velocity potential $ was shown to be 

5? + dy i+ d^ = Lagrange's equation) (1) 

where x, y, and z are the coordinates referred to axes mutually at right angles to each other. 
The velocity components in the directions of these axes are 

d<p d$ d<£ 

I assume in this paper the reader to be familiar with paper reference 1, or with tbe fundamental 
things contained therein. 

2. The configurations of velocity to be superposed for the investigation of the elementary 
technical problems of flight are of the most simple type. It will appear that it is sufficient to 
study two-dimensional flows only, in spite of the fact that all actual problems arise in three- 
dimensional space. It is therefore a happy circumstance that there is a method for the study 
of two-dimensional aerodynamic potential flows which is much more convenient for the investi- 
gation of any potential flow than the method used in reference 1 for three-dimensional flow. 
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The method is more convenient on account of the greater simplicity of the problem, there 
being one coordinate and one component of velocity less than with the three-dimensional flow. 
But the two-dimensional potential is still a function of two variables and it represents a distri- 
bution of velocity equivalent to a pair of functions of two variables. By means of introducing 
the potential a great simplification of the problem has been accomplished, reducing the number 
of functions to one. This simplification can now be carried on by also reducing the number 
of variables to one, leaving only one function of one variable to be considered. This very 
remarkable reduction is accomplished by the use of complex numbers. 

The advantage of having to do with one function of one variable only is so great, and 
moreover this function in practical cases becomes so much simpler than any of the functions 
which it represents, that it pays to get acquainted with this method even if the student has 
never occupied himself with complex numbers before. The matter is simple and can be 
explained in a few words. 

The ordinary or real numbers, x, are considered to be the special case of more general 
expressions (x+iy) in which y happens to be zero. If y is not zero, such an expression is 
called a complex number, x is its real part, iy is its imaginary part and consists of the product 
of y, any real or ordinary number, and the quantity i, which is the solution of 

■i 2 = —1; i. e.,i = V ' —1 

The complex number (x + iy) can be supposed to represent the point of the plane with the 
coordinates x and y, and that may be in this paragraph the interpretation of a complex num- 
ber. So far, the system would be a sort of vector symbolism, which indeed it is. The real 
part x is the component of a vector in the direction of the real z-axis, and the factor y of the 
imaginary part iy is the component of the vector in the y-direction. The complex numbers 
differ, however, from vector analysis by the peculiar fact that it is not necessary to learn any 
new sort of algebra or analysis for this vector system. On the contrary, all rules of calculation 
valid for ordinary numbers are also valid for complex numbers without any change whatsoever. 

The addition of two or more complex numbers is accomplished by adding the real parts 
and imaginary parts separately. 

(x + iy) + (x' +iy') = (x+x')+i(y+y') 

This amounts to the same process as the superposition of two forces or other vectors. The 
multiplication is accomplished by multiplying each part of the one factor by each part of the 
other factor and adding the products obtained. The product of two real factors is real of 
course. The product of one real factor and one imaginary factor is imaginary, as appears 
plausible. The product of ixi is taken as — 1, and hence the product of two imaginary parts 
is real again. Hence the product of two complex numbers is in general a complex number 
again 

(x + iy) (x' + iy') = (xx' - yy') +i(xy' + x'y) . 

There is now one, as I may say, trick, which the student has to know in order to get the 
advantage of the use of complex numbers. That is the introduction of polar coordinates. The 
distance of the point (x,y) from the origin (0,0) is called R and the angle between the positive 
real axis and the radius vector from the origin to the point is called <p, so that x=R cos <p; 
y=Rsm<p. Multiply now 

(i?! cos v 1 +iR 1 sin <p } ) (R 2 cos <p z +i R 2 sin ip t ). 

The result is 

R t R 2 cos <pi cos <p 2 —R 1 R 2 sin <p, sin <p 2 +i(R l R 2 cos ^ sin tp 2 +R 1 R i sin <p t cos <p 2 ) 
or, otherwise written 

RiR 2 fcosfe + <p 2 ) +i sin +<p 2 )] 
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That is: The radius R of the product is the product of the radii R t and R 2 of the two factors, 
the angle <p of the product is the sum of the angles <p t and a> t of the two factors. Further, as 
is well known, we may write 

z = R(cos <p + i sin <p) =i? e?' 
where e denotes the base of the natural logarithms. 
As a particular case 

(e = (cos <p + i sin o) n = cos n <p+i sin n <p = e rA ? 

This is Hoivre's formula. 

I proceed now to explain why these complex numbers can be used for the representation of a 
two-dimensional potential flow. This follows from the fact that a function of the complex num- 
bers, that is in general a complex number again different at each point of the plane, can be treated 
exactly like the ordinary real function of one real variable, given by the same mathematical 
expression. In particular it can be differentiated at each point and has then one definite differ- 
ential quotient, the same as the ordinary function of one variable of the same form. The 
process of differentiation of a complex function is ^determinate, in so far as the independent 
variable (x + iy) can be increased by an element (dx + idy) in very different ways, viz, in 
different directions. The differential quotient is, as ordinarily, the quotient of the increase of 
the function divided by the increase of the independent variable. One can speak of a differ- 
ential quotient at each point only if the value results the same in whatever direction of {dx+idy) 
the differential quotient is obtained. It has to be the same, in particular, when dx or dy is zero. 

The function to be differentiated may be 

F(x+iy)=R (x, y)+iT(x, y) 
where both R and T are real functions of x and y. The differentiation gives 

dx dx dx 

or again 

d£ = _.dR^dT 
idy % dy ' dy 

These two expressions must give identical results and hence are equal. That is, both the real 
parts and both the imaginary parts are equal: 

dR = _dT. bR^_dT^ 
dy dx' dy dx' 

Differentiating these equations with respect to dx and dy 

d 2 R _d 2 T_ d 2 T . d 2 T t d 2 T . 

dxdy~dy 2 ~ dx 3 ' %m dx^dy 2 

or again 

d 2 T _d 2 R_ d 2 R . d 2 R , ^R_. 0 

dxdy dy 2 dx 2 ' ' *' dx 2 ' dy 2 

Hence, it appears that the real part as well as the imaginary part of any analytical complex 
function complies with equation (1) for the potential of an aerodynamic flow, and hence can be 
such a potential. If the real part is this potential, I shall call the complex function the "poten- 
tial function " of the flow. It is not practical, however, to split the potential function in order 
to find the potential and to compute the velocity from the potential. The advantage of having 
only one variable would then be lost. It is not the potential that is used for the computation 
of the velocity, but instead of it the potential function directly. It is easy to find the velocity 
directly from the potential function. Differentiate F(x+iy) = F(z). It is seen that 

dF(z) _dR , .dT 
dz dx ' dx 
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But it was shown before that 



Hence 



bT_ 
bx~ 



bR 
'by 



dF(z) = bR_.m 
dz bz by 

The velocity has the components ^ and ^« Written as a complex vector, it would be ^ +i 
It appears therefore: 

Any analytical function F(z) can be used for the representation of a potential flow. The 

dF 

potential of this flow is the real part of this potential function, and its differential quotient 

called the "velocity function," represents the velocity at each point "turned upside . down. " 
That means that the component of the velocity in the direction of the real axis is given directly 

dF 

by the real part of the velocity function and the component of the velocity at right angles to 

dF 

the real axis is equal to the reversed imaginary part of The absolute magnitude of the ve- 

dF 

locity is equal to the absolute magnitude of -g- • 

3. I proceed now to the series of two-dimensional flows which are of chief importance for 
the solution of the aerodynamic problems in practice. They stand in relation to the straight 
line. The privileged position of the straight line rests on the fact that both the front view and 

the cross section of a monoplane are approxi- 
mately described by a straight line. The different 
types of flow to be discussed in this section have 
in common that at the two ends of a straight 
line, but nowhere else, the velocity may become 
infinite. At infinity it is zero. This suggests 
the potential function V 3 * - 1 which has discon- 
tinuities at the points ± 1 only, but it does not 
give the velocity zero at infinity. 




F=^Jz T ^l-z 

gives rise to an infinite velocity at the points 
z = ± 1 which may be regarded as the ends of the 
straight line, and in addition the velocity be- 
comes zero at infinity. A closer examination 
shows that indeed the potential function 

F=i{z-- % fz T =T) (2) 

represents the fiow.produced by the straight line extending between the points z= ±1, moving 
transversely in the direction of the negative imaginary axis with the velocity 1 in the fluid other- 
wise at rest. For its velocity function is 

F' = i- 



Fig 1.— Transverse flow, produced by a moving straight line. 



giving for points on the line a transverse velocity — 1 . This flow may be called " transverse flow." 
The velocity potential at the points of the line, i. e., for y = 0 is ^Jl—x 2 . This gives the kinetic 
energy of the flow (half the integral of the product of potential, density and normal velocity 
component, taken around the line (reference 1)). 
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r~§-*J*yj-rf<&-jr (3) 

giving an apparent mass of the straight line moving transversely equal to the mass of the fluid 
displaced by a circle over the straight line as diameter. This r emin ds us of the apparent 
additional mass of the circle itself, which is the same (ref. 1, sec. 6). It can be proved that 
the additional mass of any ellipse moving at right angles to a main axis is equal to the mass of 
the fluid displaced by a circle over this main axis as diameter (reference 6). 

The flow around the straight line just discussed can be considered as a special case of a 
series of more general flows, represented by the potential function 

F=i(z--^l) n (4) 

where n is any positive integer. n=l gives the transverse flow considered before. For n 
different from 1 the component of the transverse velocity along the straight line is no longer 
constant, but variable and given by a simple, law. Such a flow, therefore, can not be produced 
by a rigid straight line moving, but by a flexible line, being straight at the beginning and in the 
process of distorting itself. 

It is helpful to introduce as an auxiliary variable the angle 5 defined by z=cos 5. Then 
the potential function is 

F—i (cos n8—i sin n5) =i e~ inS 
where 5 is, of course, complex. The potential along the line is 

$ = sin nd (5) 
where now S is real. The velocity function is 

_,, dF dF dS n e~ inS n , . . . „. 

F =Tz=Ts' ak = -"ST— El (cos n5 ~ % sm nS) 



giving at points along the line the transverse component 
and the longitudinal component 



n sin n8 



n cos nb ,„ s 

This becomes infinite at the two ends. The kinetic energy of the flow is 

m 1 C s * sin 3 nS . „ p 

T= ~2 np J,, -^T smSdS = n *i (8) 
This impulse is given by the integral 

p y $ dx 

to be taken along both sides of the straight lines, since the velocity potential times p represents 
the impuLsive pressure necessary to create the flow. This integral becomes 



1 



ir 

sin nS sin 5 dd 



for the nth term. This is zero except for n = l. 

By the superposition of several or infinitely many of the flows of the series discussed 

F-HAfr- jtt) + A 2 (z- V*=55M-. - .+A n (s- V^3H (9) 

with arbitrary intensity, infinitely many more complicated flows around the straight line can be 
described. There is even no potential flow of the described kind around the straight line existing 
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which can not be obtained by such superposition. The kinetic energy of the now obtained by 
superposition stands in a very simple relation to the kinetic energy of the single flows which rela- 
tion by no means is self-evident. It is the sum of them. This follows from the computation of 
the kinetic energy by integrating the product of the transverse component of velocity and the 
potential along the line. This kinetic energy is 

r=|| (Ai sin 8 + A 2 sin 28 + . .) (A 1 sin 8+2A t sin 2 8 + . .)dd (10) 

But the integral 

sin nS sin m5 d8 = 0 (n^m) (11) 



is zero if m and n are different integers. For integrating two times partially gives the same 
integral again, multiplied by (mjn) 2 . In the same way it can be' proved that 



r 



cos n8 cos mS dS = 0 (n^m) (12) 



Only the squares inlntegral (10) contribute to the energy and each of them gives just the kinetic 
energy of its single term (equation (8)). 

It may happen that the distribution of the potential 3> along the line is given, and the flow 
determined by this distribution is to be expressed as the sum of flows (equation (9)). The 
condition is, for points on the line, a known function" of x is given, 

§ = A l sin 8 + A 2 sin 2d+. . . + A n sin nS + . . . (13) 

and the coefficients A are to be determined. The right-hand side of equation (13) is called a 
Fourier's series, and it is proved in the textbooks that the coefficients A can always be deter- 
mined as to conform to the condition if $ has reasonable values. At the ends 5 = 0 or t, hence <J? 
has to be zero there as then all sines are zero. 

Otherwise expressed equation (4) gives enough different types of flow to approximate by 
means of superposition any reasonable distribution of the potential over a line, with any exact- 
ness desired. This being understood, it is easy to show how the coefficients A can be found. 

Integrate ... ... 



(A t sin S + A 2 sin 28 . . .) sin n8 dS 
According to equation (11) all integrals become zero with the exception of 

sin 2 nd d8 = \ . A n 



'•X* 



Hence 

<£ sin n8 d8 (14) 



T Jo 



These values may be introduced into equation (9), and thus the potential function F is 
determined. 

Another problem of even greater practical importance is to determine the potential func- 
"tions, equation (4), which superposed give a desired distribution of the transverse component 
of velocity. The condition is now 

. sin 5 , m . sin 28 , 
u = A,-. — 1+2A,—. — r +. . . (15) 
1 sin 5 2 sin 5 

That means, now, u sin 8, a known function, is to be expanded into a Fourier's series 

u sin 8 = B X sin 8 + B t sin 28 + . . . + 5„sinn5 (16) 
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The B's may be determined by an equation like (14), and then the A's may be deduced, since 

A n =BJn (16a) 

This is always possible if the velocity component is finite along the line. These values may then 
be introduced in equation (9). 

The value of the potential function Fas given by series (13) with the values of A n substituted 
from (16) may be transformed into a definite integral which sometimes is more convenient for 
application. Let u 0 be a function of the coordinate z 0 , a point on the line joining g= — 1 and 
z = +1, and letf(z, z 0 ) be a function to be determined so that 

F=j^ f(?> Zo) • «o • dz 0 
I have found that this is satisfied by making 

f= f - [log (««-««•) -log (17) 



and this leads to a physical interpretation of u 0 . 
Hence, the velocity function 

wJl_ f +1 ¥ 

dz J-i 



dz 



The potential function and the velocity function are both thought of, then, as being the summa- 
tion of functions due to "elementary flows." An element gives rise to a potential function 

/ {z, s 0 ) u 0 dz 0 and to the velocity function dz 0 

f_df_t sin 5 a 1_ 

■ ~dz~7rshi8 cos 5 — cos 5. 



^ r 1 1 jl-i 



where the plus sign is to be taken for points on 
the positive side of the line, and the negative sign 
for those on the opposite side. In this elemen- 
tary flow, then, the velocity is parallel to the line 
at all points of the line excepting the point z 0 , 
being directed away from this point on the 
positive side of the line and toward it on the other side 




Fig. 2.— Flow around a straight line ereated by one element of the 
wing section. 



For points close to z 0 



/' u„ dz Q =^^ — 

J 9 0 T S — Z 0 

from which the value of the velocity of the flow may be deduced. If a small circle is drawn 
around the point z w it is seen that there is a flow out from the point 2 0 of amount u 0 dz a per sec- 
ond on the positive side and an inflow of an equal amount on the other side; so that this is 
equivalent to there being a transverse velocity u 0 at points along the element dz m positive on 
one side, negative on the other. The total flow around the line due to f{z , z 0 ) u 0 dz 0 is illustrated 
in Figure 2. 

Substituting the value off in F' 



F , = , r+ i i_ j 1-2 
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Therefore, for any point on the real axis, the transverse Telocity is u 0 and the longitudinal 
velocity 



■2 3 



IT Z — Z 0 V 1 —2 

Or, interchanging symbols, writing z for z 0 and vice versa 



f +1 I wfe /J-s 2 



(19) 



For a point near the edge on the positive side, write z 0 = l — e 



or, substituting -yj2e=sm 5 edl 



Tiffin J^^Vrii 



-2 



(19a) 



For the discussion of the elements of the wing theory, in addition to the flows mentioned, 
there is one flow which needs a discussion of its own. This is given by the potential function 



The velocity function of this flow is 



F=A 0 sin -1 z 



(20) 



'■yJT=Z>~ 

I shall call this flow "circulation flow" as it 
represents a circulation of the air around the 
line. The transverse component of the velocity 
at points along the line is identically zero. 

The circulation flow does not quite fit- in 
with the other ones represented by equation 
(4), because the potential function (20) is a 
multiple valued one, the values at any one 
point differing by 2tt or multiples thereof. All 
this indicates that the flow is a potential flow, 
it is true, but it does not conform to the condi- 
tion of a potential flow when considered as in 
equation (4). 

This is in accordance with the physical 
consideration, that it is impossible to produce 
this flow by an impulsive pressure over the 
straight line. Such a pressure would not per- 
form any mechanical work, as the transverse 
components of velocity at points along the line are zero. Tbe kinetic energy of this flow, on 
the other hand, is infinite, and hence this flow can not even be completely realized. Still it 
plays the most important part in aerodynamics. 

The best way to understand this flow and its physical meaning is to suppose the line to be 
elongated at one end, out to infinity. On the one side the potential may be considered zero. 
Then it is constant and will be equal to 2x on the other side. The transverse velocity component 
is finite. Hence the flow can be produced by a constant impulsive pressure difference along 
this line extending from the edge 2 = 1 to infinity. This pressure difference makes tbe fluid 
circulate around the original straight line, the pressure along tbe line itself being given by the 
potential function (20) and not performing any work. 




Fig. 3.— Circulation flow around a straight line. 
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A pressure difference along an infinite line does never actually occur. At least t does not 
occur simultaneously along the whole line. A very similar thing, however, occurs very often 
which has the same effect. That is a constant momentum being transferred to the air at 
right angles to an infinite straight line at one point only, but the point traveling along the line, 
so that the final effect is the same as if it had occurred simultaneously. This is the fundamental 
case of an airplane flying along that infinitely long line. During the unit of time it may cover 
the length V and transfer to the air the momentum L, equal to the lift of the airplane. Then 
the impulse of the force, per unit of length of the line, is L/V, and hence the potential difference is 

y^- That makes A 0 in equation (20) A 0 = ^y^ - If the airplane has traveled long enough, 

the flow in the neighborhood of the wing, or rather one term of the flow, is described by the 
circulation flow, provided that the airplane is two-dimensional, that is, has an infinite span. 
The velocity at the end of the wing 2= due to this circulation flow 

F=A a sin -1 2, 

where 
is 

F ^ e = V^A=o = &rpF (sin 5) 5=0 (21a) 
II. THEORY OF WING SECTION. 

4. The investigation of the air flow around wings is of great practical importance in view 
of the predominance of heavier-than-air craft. It is necessary to divide this problem into two 
parts, the consideration of the cross section of one or several wings in a two-dimensional flow, 
and the investigation of the remaining effect. This chapter is devoted to the first question. 

All wings in practice have a more or less rounded leading edge, a sharp trailing edge and 
the section is rather elongated, being as first approximation described by a straight line. The 
application of the aerodynamic flows around a straight line for the investigation of the flow 
around a wing section suggests itself. I have shown in section (3) how the potential flow 
around a straight line is determined, for instance, from the transverse components of velocity 
along this line. Only one type of flow, the circulation flow, is excepted. This flow does not 
possess any transverse components at the points of the line and hence can be superposed on a 
potential flow of any magnitude without interfering with the condition of transverse velocity. 
I have shown, on the other hand, that it is just this circulation flow not determined so far, 
which gives rise to the chief quantity, the lift. It is, therefore, necessary to find some additional 
method for determining the magnitude of the circulation flow. 

This magnitude of the circulation flow is physically determined by the fact that the air 
is viscous, no matter how slightly viscous it is. The additional condition governing the magni- 
tude of the circulation flow can be expressed without any reference to the viscosity and was 
done so in a very simple way by Kutta. The condition is very plausible, too. Kutta's condition 
simply states that the air does not flow with infinite velocity at the sharp, rear edge of 
the wing section. On the contrary, the circulation flow assumes such strength that the air 
leaves the section exactly at its rear edge flowing there along the section parallel to its mean 
direction. The wing as it were acts as a device forcing the air to leave the wing flowing in a 
particular direction. 

Consider, for instance, the wing section which consists merely of a straight line of the length 2 . 
The angle of attack may be a. The flow produced by this line moving with the velocity V is 
then represented by the potential function 

F= V sin a • i e~ si 

which gives a constant transverse component of velocity along the wing, as shown in equation 
(6) for n = l. The real axis is parallel to the straight line, its origin is at the center of the line. 
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The infinite longitudinal velocity at the rear end is 

— Fsin 



(sin 5)3=0 

The angle of attack may now be assumed to be small and I change slightly the way of 
representing the flow, turning the real axis of coordinates into the direction of motion. Instead 
of referring the flow to the line really representing the wing section I consider the straight line 
between z— ±1, which differs only slightly from the wing and is parallel to the motion. The 
transverse components of the flow relative to this line are approximately equal to the transverse 
velocity relative to. the wing section at the nearest point and therefore constant again and equal 
to Fsin a. Therefore, this way of proceeding leads to the same flow as the more exact way 
before. It also gives the same infinite velocity at the rear end. 

This velocity determines the magnitude of the circulation flow 

F=i 0 sin- 1 2 (21) 

by the condition that the sum of their infinite velocities at the edge is zero. 

V sin a , / —A 0 , } = 0 
(sin 5)j=o 0 (sm 5)s=o 

and hence A a = V sin a. The lift is therefore 

£=#tt V 2 p sin a. 

The lift coefficient, defined by Ol = — — — , since the chord = is therefore, 

Oi = 2t sin a, or approximately 2tra (22) 

and 

L=7 s |S£xa (23) 

where S denotes the area of the wing. 

The representation of the flow just employed is approximately correct and gives the same 
result as the exact method. This new method now can be generalized so that the lift of any wing 
section, other than a straight line, can be computed in the same way, too. The section can be 
replaced with respect to the aerodynamic effect by a mean curve, situated in the middle between 
the upper and lower curves of the section and having at all points the same mean direction as 
the portion of the wing section represented by it. The ordinates of this mean wing curve may 

be £, the abscissa z, so that the direction of the curve at each point is This direction can 

be considered as the local angle of attack of the wing, identifying the sine and tangent of the 
angle, with the angle itself. Accordingly it is variable along the section. Since the velocity of 
the air relative to the wing is approximately equal to the velocity of flight, the component at 

VdP 

right angles to the s-axis is As before, the infinite velocity at the rear edge is to be found. 

It is, according to equation (19 a) 

7T (sm §){=„ J-i dz -\/ / -x ' 
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At the rear edge x = l. The mean apparent angle of attack, that is the angle of attack of the 
straight line giving the same lift as the wing section, is found by the condition that this infinite 

"V sin & i 

value must be the same as that deduced for a straight line; viz, ■ = — -z- * Hence, replacing 

sin o 

sin a by a 

a ^J_C^^f^Bdx length = 2 (25) 
a" J- i dx\ 1-x 

«--§/:: ivsi^- 1 

This formula holds true for any small angle of attack of the section. The integral can now 
be transformed into one containing the coordinate | rather than the inclination ^ of the wing 
curve, provided that the trailing edge is situated at the ar-axis, that is, if | is zero at the end 
x= +1. This transformation is performed by partial integration, considering g&asa factor 
to be integrated. It results 

«'=- P 1 u W/ * length = 2 (26) 
tJ-, V-x) -\I1-x 2 



4 &e 



(26a) 



The important formula (26) gives the mean apparent angle of attack directly from the 
coordinates of the shape of the wing section. The mean height % of the section has to be inte- 
grated along the chord after having been multiplied by a function of the distance from the 
leading edge, the same for all wing sections. This integration can always be performed, whether 
the section be given by an analytical expression, graphically or by a table of the coordinates. 
In the latter case a numerical integration can be performed by means of Table II, taken from 
reference 4. The figures in the first col umn give the distance from the leading edge in per cent 
of the chord. The second column of figures gives factors for each of these positions. The 
height £ of the mean curve of the section over its chord, measured in unit of the chord terms, 
is to be multiplied by the factors, and all products so obtained are to be added. The sum gives 
the apparent angle of attack in degrees. 

5. The lift of a wing section as computed in section (4) is caused by the circulation flow 
s ymm etrical with respect to the straight line representing the wing. Hence the pressure creat- 
ing this lift is located symmetrically to the wing, its center of pressure is at 50 per cent of the 
chord, it produces no moment with respect to the middle of the wing. This lift is the entire 
lift produced by the wing. It is not, however, the entire resultant air force. The remaining 
aerodynamic flow in general exerts a resultant moment (couple of forces) and this moment 
removes the center of pressure from its position at 50 per cent. 

If the wing section is a straight line of the length 2, its apparent transverse mass is xp, as 
seen in section (4). The longitudinal mass is zero. Hence, according to reference 1, the result- 
ant moment is 

M= V 2 | x sin 2a length =2 (27) 

M~ V*^2xa length = 2 (28) 

Both the exact and the approximate expressions give the constant center of pressure 25 per cent 
of the chord from the leading edge, as results by dividing the moment by the lift (23). 

The straight sections considered have a constant center of pressure, independent of the angle 
of attack. The center of pressure does not travel. This is approximately true also for symmet- 
rical sections with equal upper and lower curves, where the center of pressure is also at 25 
per cent. If, however, the upper and lower curves are different and hence the mean section 
52201—25 11 
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curve is no longer a straight line, the potential flow produced at the angle of attack zero of the 
chord not only gives rise to the circulation flow and thus indirectly to a lift, but also creates a 
moment of its own. It is simple to compute this moment from the potential flow, which is 
represented in equation (9) as a superposition of the flows, equation (4). 

The longitudinal velocity relative to the line is, according to equation (7), 

/ , cos 5 , . 2 cos 28 , , . n cos nS\ , Tr 

As the section is supposed to be only slightly curved, ^ is always small, so are, therefore, the 

coefficients A n when compared to V, so that they may be neglected when added to it. The 
pressure at each point along the line, according to reference 1, is 

t-i* 

The present object is the computation of the resultant moment. When really forming the 
square of the bracket in the last expression, the term with Vindicates a constant pressure and 
does not give any resultant moment. The squares of the other terms are too small and can be 
neglected. There remains only the pressure, 

p-pFfj, cogs 3 cos 35 \ 

c \ 1 sin S 2 sin § / 

giving the resultant moment about the origin 

*r a rr f V . COS 5 , . # COS 25 \ . . , 

M=2 V P J f (A zTs+A, ) cos 8 sin ids, 

since the density of lift is twice the density of pressure, the pressure being equal and opposite 
on both sides of the wing. But according to (12) 



cos nd cos m8 d8= 0 (12) 
if m and n are different integers. Hence there remains only one term. The resultant moment is 

M=2 P VA,f 

A t was found according to equation (14) by means of the integral 



dx 

Hence the moment is 



i" 

*■ Jo 



V$-sitf8d5 



dx 



sin 2 SdS 



or, expressed by x 

M=2pV*f"f x TjT=x~*dx (29) 

By the same method as used with integral (25) this integral can be transformed into 

M=B P V* ( +1 4== (30) 

J-l -yfi—X 2 



ELEMENTS OP WING SECTION THEORY AMD WISG THEOBY. 



153 



It has been shown that for a chord of length 2, the center of pressure has a lever arm ^ and the 

lift is F 2 § • 2t<x ■ 2, giving a moment F 2 pxa; so that an angle of attack corresponds to a 

moment F 2 pira. Consequently the resultant moment is the same as if the angle of attack is 
increased by the angle 

""'U-!&^ tb=2 <31) 



1/2 —43? 



length =1 (32) 



It is readily seen that this angle is zero for sections with section curves equal in front and 
in rear. Hence such sections have the center of pressure 50 per cent at the angle of attack zero 
of the mean curve, that is, for the lift (24) produced by the shape of the section only. The 
additional lift produced at any other angle of attack of the chord and equal to the lift as produced 
by the straight line at that angle of attack has the center of pressure at 25 per cent. Hence a 
travel of the center of pressure takes place toward the leading edge when the angle of attack is 
increased, approaching the point 25 per cent without ever reaching it. The same thing happens 
for other sections with the usual shape. At the angle of attack zero of the chord the lift pro- 
duced was seen to be 2xY^pa ', i. e., from (26) 

r - 172 £ / f +1 dx £ 

and the moment, see equation (30), 

M=2pV* ( +1 ^Jt_ (30) 
j-i V-i — 

giving the center of pressure at the distance from the middle 



x%dx 



—I {.l-X)-yJT=X* 

The lift produced by the angle of attack of the chord, equation (23) as before has the center of 
pressure 25 per cent. The travel of the center of pressure can easily be obtained from this 
statement. The moment about the point 25 per cent is independent of the angle of attack. 
The center of pressure in ordinary notation at the angle of attack zero is 

" +1 x%dx 



r 

GP = 50%-50 J- 1 % 
J-i (J-a;)yi-a? 



The computation of the mean apparent angle of attack with respect to the moment is done in 
the same way as that of the angle with respect to the lift. Table II, gives the coefficients for 
numerical integration, by means of two ordinates only, to be used as the other figures in Table II. 
The final sum is the mean apparent angle in degrees. 

6. The problem of two or even more wing sections, combined t o a biplane or multiplane and 
surrounded by a two-dimensional flow can be treated in the same way as the single wing section. 
The two sections deter mine by their slope at each point a distribution of transverse velocity 
components along parallel lines. The distribution determines a potential flow with a resultant 
moment. According to Kutta's condition of finite velocity near the two rear edges, the potential 
flow in its turn determines a circulation flow giving rise to a lift and moment. The physical 
aspect of the question offers nothing new, it is a purely mathematical problem. 
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This mathematical problem has not yet been solved in this extension. I have attacked 
the problem within a more narrow scope (reference 4) . The method followed by me amounts to 
the following considerations: 

Equation (13) represents different types of Sow around one straight line, consisting in a 
motion of the air in the vicinity of the straight line only. Now the motion of the flows with high 
order n is more concentrated in the immediate neighborhood of the straight line than the flows 
of low order n. The transverse velocity components along the line, determining the flow, change 
their sign (n-1) times along the line. With large n, positive and negative components follow each 
other in succession very rapidly so that their effect is neutralized even at a moderate distance. 

Hence the types of flow of high order n around each of a pair of lines will practically be the 
same as if each line is single. The flows of high order do not interfere with those of a second line 
in the vicinity even if the distance of this second line is only moderate. It will chiefly be the 
types of flow of low order, the circidation flow n=0, the transverse flow n = l, or it may be the 
next type n = 2 which differ distinctly whether the wing is single or in the vicinity of a second 
wing. Accordingly, I computed only the flows of the order n = 0 and n = l, the circulation flow 
and the transverse flow for the biplane and used the other flows as found for the single section. 

The results are particularly interesting for biplanes with equal and parallel wings without 
stagger. Their lift is always diminished when compared with the sum of the lifts produced by 
the two wings when single. The interference is not always the same. If the sum of the angle 
of attack and the mean apparent angle of attack with respect to the moment is zero, or other- 
wise expressed, at the angle of attack where the center of pressure is at 50 per cent, it is par- 
ticularly small. The lift produced at the angle of attack zero is diminished only about half as 
much as the remaining part of the lift produced by an increase of this angle of attack. 

This second part of the lift does not have its point of application exactly at 25 per cent of 
the chord, although its center of pressure is constant, too. This latter is quite generally valid for 
any two-dimensional flow. At any angle of attack zero arbitrarily chosen, the configuration of 
wing sections produces a certain lift acting at a certain center. The increase of the angle of 
attack produces another lift-acting at another fixed point. Hence the moment around this 
second center of pressure does not depend on the angle of attack; and the center of pressure at 
any angle of attack can easily be computed if the two centers of pressures and the two parts of 
the lift are known. 

The resultant moment of the unstaggered biplane consisting of portions of equal and 
parallel straight lines is again proportional to the apparent transverse mass, as the longitudinal 
mass is zero (reference 1). This mass is of use for the considerations of the next chapter, too. 
Therefore, I wish to make some remarks concerning its magnitude. If the two straight lines are 
very close together, the flow around them is the same as around a line of finite thickness and is 
almost the same as around one straight line. Its apparent mass is the same, too, but in addition 
there is the mass of the air inclosed in the space between the two lines and practically moving 
with them. Hence the mass is approximately 

b{bj + Ji) P 

where 5 is the length of the lines and h their distance apart, if the distance Ti of the lines is small. 
For great distance, on the other hand, the flow around each of the lines is undisturbed, the 
apparent mass is twice that of the flow around each line if single. It is therefore 

2V-.P 
4 

For intermediate cases the apparent mass has to be computed. Particulars on this computa- 
tion are given in reference 4. Table 1 gives the ratio of the apparent mass of a pair of lines to 

that of one single line for different values of |« This ratio, of course, is always between 1 and 2. 
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These few remarks on the theory of biplane sections seem to be sufficient in this treatise 
on the elements of wing theory. The student will find full information on the subject in my 
paper on biplanes, reference 4. The remarks laid down here, I hope, will assist him in under- 
standing the leading principles of the method there employed. 

HI. AERODYNAMIC INDUCTION. 

7. The last chapter does not give correct information on the aerodynamic wing forces, 
since the flow in vertical longitudinal planes was supposed to be two-dimensional. The, vertical 
layers of air parallel to the motion were supposed to remain plane and parallel and only the 
distortion of the two other planes at right angles to it was investigated. This is a very incom- 
plete and arbitrary proceeding, for the vertical longitudinal layers do not remain plane, as 
little as any other layers remain plane. It is therefore necessary to complete the investigation 
and to assume now another set of layers, parallel to the lift, to remain plane, thus studying 
the distortion of the vertical longitudinal layers. Accordingly. I will now assume that all 
vertical layers of air at right angles to the motion remain plane and parallel, so that the air 
only moves at right angles to the direction of flight. Hence, I have now to consider two- 
dimensional transverse vertical flows. This consideration, it will appear, gives sufficient 
information on the motion of the air at large, whereas the preceding investigation gives infor- 
mation on the conditions of flow in the vicinity of the wing. Both, the longitudinal two- 
dimensional flow studied before and the two-dimensional flow to be studied presently, possess 
vertical components of velocity. Both flows and in particular these vertical components 
are to be superposed, and thus one can determine the final aerodynamic pressures and resultant 
forces. 

The transverse vertical layer of air is at rest originally. The wings, first approaching it, 
then passing through it and at last leaving it behind them, gradually build up a two-dimensional 
flow in each layer. The distribution of impulse creating this flow is identical with the distri- 
bution of the lift over the longitudinal projection of the wings. It is immaterial for the final 
effect whether all portions of the wings at every moment have transferred the same fraction 
of the momentum to a particular layer or not. The final effect and hence the average effect 
is the same as if they always have. They actually have if all wings are arranged in one trans- 
verse plane — that is, if the airplane is not staggered. It may be assumed at present that at 
each moment each layer has received the same fraction of the impulse from every portion of 
the wings and it follows then that the shape of the co nf iguration of the two-dimensional flow 
is always the same and that it is built up gradually by in-creasing its magnitude while not 
changing its shape, beginning with the magnitude zero at a great distance in front of the wing 
and having obtained its final magnitude at a great distance behind the wings. 

The potential of the final two-dimensional flow long after the wings have passed through 
the layer is easy to find, for the impulsive pressure creating it is known along the longitudinal 
projection of the wings. It is identical with the distribution of the lift over this projection, 
acting as long as the airplane stays in the layer. This is the unit of time, if the thickness of 
the layer is equal to the velocity of flight. Hence the potential difference along the longi- 
tudinal projection of the wings is equal to the density of the lift along this projection divided 
by the product of the density of air and the velocity of flight, since the velocity potential is 
equal to the impulse of the pressure creating the flow, divided by the density. In general the 
longitudinal projections of the wings can be considered as lines. The density of lift per unit 
length of these lines is then equal to the difference of pressure on both sides, and hence the 
density of the lift is proportional to the difference of the potential on both sides. This state- 
ment determines completely the final two-dimensional flow in the transverse vertical layer, 
and nothing remains unknown if the distribution of the lift over the wings is given. The 
actual determination of the flow is then a purely mathematical process. 
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For the present purpose, however, not the final transverse flow but the vertical flow at the 
moment of the passage of the wings is of interest. It is this flow that is to be superposed on 
the longitudinal flow in order to determine the actual air forces. It has already been men- 
tioned that this flow can be supposed to differ from the final flow in magnitude only. It 
remains therefore only to find the ratio of momentum already transferred while the wing passes 
through the layer, to the momentum finally to be imparted. 

The fraction X A seems to me more plausible than any other fraction. The effect of the 
wing on the layer is the same at equal distances from the layer, whether in front or back of it 
and this would involve the factor }4- K~is not necessary, however, to have recourse to a mere 
assumption in this question, however plausible it may be. It can be proved that the assump- 
tion of %, is the only one which does not lead to a contradiction with the general principles of 
mechanics. I proceed at once to demonstrate this. 

If the transverse flow in the plane of the wings is found, only the vertical component down- 
ward u', called the induced downwash, is used for the computation. This downwash can be 
positive or negative, but in general is positive. Such downwash in the neighborhood of a 
portion of wing changes the motion of the air surrounding the wing portion relative to it. The 
induced downwash is always small when compared with the velocity of flight. Hence, its 
superposition on the velocity of flight at right angles to it does not materially change the mag- 
nitude of the relative motion between the wing and the air in its vicinity. It changes, however, 
the direction of this relative velocity, which is no longer parallel to the path of the wing but 
inclined toward the path by the angle whose tangent is u'JV. This has two important 
consequences. 

The flow produced and hence the air force no longer correspond to the angle of attack be- 
tween the wing and the path of flight but to the angle given by the motion of the wing relative 
to the surrounding portion of the air. In most cases the angle of attack is decreased and the 
effective angle of attack, smaller than the geometric angle of attack between path and wing, 
determines now the flow and the air forces. Hence, the lift" in general is smaller than would 
be expected from the geometric angle of attack. The angle of attack in the preceding chapter 
on the wing section is not identical with the geometric angle between the chord and the direc- 
tion of flight but with the effective angle of. attack, smaller in general, as there is an induced 
downwash motion in the vicinity of the wing. Therefore the geometrical angle of attack is 
decreased by 

<H=y (34) 

That is not all. The lift is not only decreased but its direction is changed, too. It is no 
longer at right angles to the path of flight, but to the relative motion between wing and adjacent 
portion of air. It is generally turned backward through an angle equal to the induced angle 
of attack. The turning backward of the lift by itself does not materially change the magni- 
tude of the lift, as the angle is always small; the vertical component of the lift remains almost 
the same, but the effective angle of attack has to be decreased. In addition to this the air force 
has now a component in the direction of the motion. The wing experiences an "induced" 
drag, in addition to the drag caused by the viscosity of the air, not discussed in this paper, 
and the induced drag is often much larger than the viscous drag. The density of the induced 

drag is cLL y where dL is the density of lift, as can be directly seen from Figure 4. 

dD t = dL^ (35) 

The existence of a drag could have been anticipated, as there must be a source of energy 
for the creation of the transverse flow under consideration. The final kinetic energy of this 
flow in a layer of thickness V is 

J ~dL u 
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and this energy is to be delivered by the wing per unit of time, as during this unit of time 
another layer has been put into motion in the way discussed. On the other hand the energy 
delivered by the wings is the integral over the drag multiplied by the velocity, that is, f&L u'. 

From which follows immediately u' =^ u, and it is thus confirmed that the transverse flow 

is only half formed when the wings are passing through the vertical layer. 

8. The problem is thus solved in general if the shape of the wings and the distribution 
of lift over the wings is known. Before passing to special wing arrangements and distributions 
of lift, in particular to the simple monoplane, there is one general problem to be discussed. 
The longitudinal projection of the wings being given, as well as the entire lift, the induced 
drag depends on the distribution of the lift over the projection. The drag is desired to be as 
small as possible. The question arises, What is the distribution of lift giving the smallest 
induced dragl The importance of this question is at once obvious. 

The entire lift and the entire induced drag of the wings are found again as important 
characteristics of the final transverse flow, discussed in the last section. The resultant lift 
is equal to the resultant vertical momentum of this flow for the thickness of the layer equal 
to the velocity V, and the induced drag is equal to the kinetic energy in the same layer divided 
by V. The problem is therefore to find such a two-dimensional flow produced by impulsive 
pressure over the longitudinal projection of the wings as possesses a given magnitude of the 
vertical momentum, and the kinetic energy of which i3 a minimum. 

It is sufficient for elementary questions to consider only arrangements of wing symmetrical 
with respect to a vertical longitudinal plane, giving moreover horizontal lines in the longitudinal 
projections. The results are valid for all conditions (reference 1). It is then easy to find the 
solution. The momentum of several flows superposed on each 
other is the sum of their single momenta. The flow is of the de- 
sired kind if the superposition of any otherflow with the resultant 
vertical momentum zero increases the kinetic energy of the flow. 

The velocity of the superposed flow can be assumed to be 
small, for instance, so that its own kinetic energy, containing the 
square of the velocity, can be neglected. The impulsive pressure 
along the projection of the wings necessary to create the super- 
posed flow acts along a path deter min ed by the magnitude of the 

downwash at the same points. The increase of kinetic energy is — ~ ~~ *| u' 

~p P S u $ (fob fiG. i— Diagram showing the creation at the 

* induced drag. 

where f # dz = 0. 

It is readily seen that the first expression can be identically zero for any distribution of the 
potential <f> restricted by the second condition only if the downwash u' is constant over the 
entire projection of the wings. Only then a transfer of a portion of lif t from one point to another 
with smaller downwash is impossible, whereas this proceeding in all other cases would lead to a 
diminution of the induced drag. It is thus demonstrated: 

The induced drag is a m in i m um, if the transverse two-dimensional flow has a constant 
vertical component of velocity along the entire projection of the wings. 

For wings .^without stagger it follows then that the induced angle of attack u'/V is 
constant over all wings. 

The magnitude of the minimum induced drag of a system of wings is easily found from the 
apparent mass pK of their longitudinal projection in the two-dimensional transverse flow. 
For the vertical momentum equal to the lift is uV pK=L where u is the constant downwash 
of the final flow. This gives 

L 



L so* 



■90° \ Direction of wind 
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The induced drag is equal to the kinetic energy divided by V 



It follows therefore that the minimum induced drag is 




(36) 




(37) 



J5T is a constant area determined by the longitudinal projection of the wings. It is the area of 
the air in the two-dimensional flow having a mass equal to the apparent mass of the projection 
of the wings. 

The results (36) and (37) show that the minimum induced drag can be obtained from the 
consideration that the lift is produced by constantly accelerating a certain mass of air downward 
from the state of rest. The apparent mass accelerated downward is at best equal to the apparent 
mass of the longitudinal projection of the airplane in a layer of air passed by the airplane in the 
unit of time. 

In practical applications the actual induced drag can be supposed to be equal to the mini- 
mum induced drag, and the average induced angle of attack equal to (37). It is, of course, 
slightly differenty but the difference is not- great as can be expected since no function changes 
its value much in the neighborhood of its minimum. 

I proceed now to the application of the general theory of induction to the case of the mono- 
plane without dihedral angle, giving in longitudinal projection a straight line of the length b. 
Consider first the distribution of lift for the minimum induced drag. It is characterized by the 
transverse potential flow with constant vertical velocity component along this straight line. 
This flow has repeatedly occurred in the earlier parts of this paper. For the length 2 of the line, 
it is the transverse flow given by equation (2) or by equation (4) and n = l. The potential 
function for the length h is 



The density of the lift per unit length of the span is equal to the potential difference of the final 
flow on both sides of the line, multiplied by Fp. 



where cos 5=-y Plotted against the span, the density of lift per unit length of the span is 

represented by half an ellipse, the multiple of sin 5 being plotted against cos 5. The lift there- 
fore is said to be elliptically distributed. 




(38) 



giving the constant vertical velocity component along the line 





(39) 




Hence 
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The apparent mass of the line with the length b according to equation (3) is equal to 

4 

Hence, with this distribution of lift, the minimum induced drag is, according to equation (36) 

ytv-2 p 

and the constant induced angle of attack according to equation (37) is 

L (41) 



The density of lift (39) per unit length of the span together with the chord c, different in general 
along the span, and with equation (22) determines the effective angle of attack at each point, 
including the apparent mean angle of attack of the section. Namely, from equation (39), 

dL 

^ _ density of lift dx _ ££ sin 5 ^2) 
2t chord Y- 1 SkcY* | Sir'cPg 

The geometric angle of attack is greater by the constant induced angle of attack, and hence 



_ 2L sin S , L f * , ire 



-«( 1+ W3srd (43) 



Equation (40) indicates the importance of a sufficiently large span in order to obtain a small 
induced drag. 

Any distribution of lift ^ over the span other than the elliptical distribution is less simple 

to investigate, as then the induced downwash is variable. The distribution of lift gives directly 
the distribution of the potential difference along the two-dimensional wing projection. 

dL 

A$ = ^ (44) 
Vp 

The transverse two-dimensional flow can now be obtained by superposition of types of flow 

2x 

given by equation (4) with z = y > as now the length of the line is not 2 but 6. The condition 
that the superposition of such flows gives the required potential difference, viz, 

dL 



1 dx 

g A§ = jy-p =A sin 5+ A 2 sin 25+ • • - +A n sin nS+ (45) 

Hence, the distribution of the density of lift, divided by 2 Vp is to be expanded into a Fourier's 
series. The induced angle of attack results then, according to equation (15), 

or { =^y|g-^ (j4, sin S-r2A 2 sin 25-5- • • -fra A n sin nS + , ■ ■ •) (46) 



